Generation of spatiotemporal correlated noise in 1+1 dimensions 
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We propose a generalization of the Ornstein-Uhlenbeck process in 1 + 1 dimensions which is 
the product of a temporal Ornstein-Uhlenbeck process with a spatial one and has exponentially 
decaying autocorrelation. The generalized Langevin equation of the process, the corresponding 
Fokker-Planck equation, and a discrete integral algorithm for numerical simulation is given. The 
process is an alternative to a recently proposed spatiotemporal correlated model process [J. Garci'a- 
Ojalvo et al., Phys. Rev. A 46, 4670 (1992)] for which we calculate explicitely the hitherto not 
known autocorrelation function in real space. 
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I. INTRODUCTION 

Noise induced phenomena are subject of considerable 
recent attention 0, 01 ■ After considering in the early 
phase systems with only few degrees of freedom in the 
last decade the effects of noise in spatially distributed sys- 
tems have been investigated ^ . In this context stochas- 
tic model processes are necessary to mimick spatiotempo- 
ral fluctuations of different origin. If characteristic time 
and length scales of system and noise are clearly sepa- 
rated, the use of a spatiotemporal Gaussian white noise 
may be justified, but it can also lead to spurious results, 
1^ |j> 1^ ' There are physical situations where the char- 
acteristic scales are not well separated, e.g. in externally 
driven systems Jij, or where the square of the driving 
stochastic process is involved. Both is the case in elec- 
trohydrodynamic convection in nematic liquid crystals, 
cf. e.g. I3, driven by external stochastic electric fields 



Further examples are the influence of spatiotempo- 
ral colored noise on spatiotemporal chaos modeled by the 
complex Ginzburg-Landau equation 10] and on networks 
of excitable systems displaying spatiotemporal stochas- 
tic resonance llj. This sufficiently motivates to model 
correlated spatiotemporal fluctuations. For an approach 
based on different grounds see, e.g., 0. 

A frequently used spatiotemporal correlated model 
process was introduced by Garcia-Ojalvo et al. (GSR) 
[13I ] who considered in spatial dimension d, r G i?'' the 
stochastic partial differential equation (PDE) 



r^^(r,t) 



(l-A2A)(p(r,t)-be(r,i), (1) 



where the additive driving process ^,{^^t^ is Gaussian 
distributed with zero mean and with autocorrelation 
i^«(r, t- r', t') = (e(r, t)e(r', t')) = a''8{v ~ r')6{t - f). 

The heuristics of Eq. is evident: the diffusive term 
effectively reduces the life-time of Fourier components 
with wavelengths short compared to A, see also below. 
For A = it reduces to the Langevin equation defin- 
ing the common temporal Ornstein-Uhlenbeck process 
(OUP), 0,111, E3- The solutions are thus, in a sense, 
generalizations of the OUP. Equation ^ belongs to a 



class of stochastic PDEs for which existence and unique- 
ness of the solution are proven rigorously [l7l ]. It is dis- 
cussed also in the context of reaction-diffusion systems, 
see e.g. |T^ . and within a generating functional approach 

In this paper we propose an alternative spatiotempo- 
ral generalization of the OUP in 1+1 dimension which 
is simply the product of a temporal OUP with a spatial 
one and has exponentially decaying autocorrelation. To 
make the paper self-contained and to introduce the no- 
tation which is used in the sequel we shortly recall in 
Sec. II basic properties and numerical generation of the 
common OUP in one (temporal) dimension. The scaling 
necessary when transforming between discrete and con- 
tinuous formulations is carefully discussed. In Sec. Ill the 
generalized OUP is constructed independently within a 
spatially discretized scheme and in a continuous version 
as the solution of a stochastic PDE different from Eq. 
(0. Subsequently, conditions which ensure stationarity 
and homogeneity are discussed, the generalized Fokker- 
Planck equation and its stationary solution are given, and 
numerically generated data are compared with the ana- 
lytic results. In Sec. IV the autocorrelation function of 
the GSR process is explicitely calculated in real space for 
d = 1 in both continuous and discrete formulation and 
compared with numerical results. Previous work studied 
the behaviour in real space only for spatial dimensions 
d = 2, [Hill], and d = 3, ^19], cf. however 20]. Con- 
trary to the folklore |10l 122 ] , the autocorrelations of the 
GSR process decay not exponentially but in a more in- 
tricate way. 

Problems connected with the generalization to higher di- 
mensions are shortly discussed in the concluding section. 



II. THE ORNSTEIN-UHLENBECK PROCESS 



The OUP ist the only stationary Gaussian Markov pro- 
cess with exponentially decaying autocorrelation (Doob's 
Theorem [23 )• Realizations ri{t) of the OUP can be gen- 



erated solving the Langevin equation 



(2) 



where ^(t) is a Gaussian white noise with {£,{t)(,{t')) = 
a^6{t — t'). In mathematicany precise form Eq. ||2Jl reads 



Tdr]{t) = -T]{t)dt + dW{t), 



(3) 



where W is a Wiener process with {W{t)W{t')) = 
al min(t,t'); note dW{t)/dt = £,{t). Solving Eq. ® with 
initial condition ?](to) = Vo gives 



■q{t) = 77oe 



dW(s)e'''\ (4) 



which has the autocorrelation 

^2 ^ 



^ ^ (5) 

The process becomes stationary if the initial values are 
Gaussian distributed with zero mean and variance cr^ /2r, 
or in the limit i,t' — > c», or for ~* — oo; it is then the 
OUP. We denote the stationary part of the autocorrela- 
tion function as 



K'Ht-t') ^ -^e 



-\t~t'\/r 



(6) 



Naturally, K'^{t — t') solves the inhomogeneous equation 
which is obtained by multiplying Eq. Q with ri{t) given 
by Q and averaging 

r^K^{t-t') = -if''(<-i') + e(i'-t)— e^(*'-*)/". (7) 
dt T 



1 pt+At 

7]{t + At) = r]{t)e-^*'^ + _e-(*+A*)/- / dW{s)e''\ 



Given the value ri{t) we can obtain + At) as 

/.t+At 
,-(t+At)/r I 

(8) 

where the last term of the r.h.s is a stochastic increment. 
The increments in non-overlapping time intervals are ob- 
viously independent; they have zero mean and the vari- 
ance, cf. e.g. 0, 



^ f 1 - 
2t V 



(9) 



Introducing the notation rjt — ri{At ■ t), where At is fixed 
and t — . . . , —1, 0,1,.. ., one obtains, for every choice of 
At, an exact recursion relation for equidistant discrete 
times (discrete integral algorithm), 



rit+i = rjt e 



/2t 



t+i, 



(10) 



where are independent zero mean Gaussian random 
numbers with variance one, cf. |25| . 



For small At a Taylor expansion of the r.h.s. of Eq. ^ 
leads to 

r/(i + At) - 7j{t) (^1-^^+ ^e(i + At). (11) 

The discrete version of Eq. (|lll) is obtained using the 
above notation and replacing ^(t) by crt^t/\/A<. This en- 
sures the correct autocorrelation in the continuum limit 
observing limAt^o ^t,t' /At — 6(t — t') and amounts to a 
rescaled variance — a'^/At, 



Vt+i = ?7t 1 



At. 



-o-t^+i- 



(12) 



All results are consistent: Eq. Hll() can be derived from 
Eq. (O using an Euler discretization, and Eq. H12(l from 
Eq. HlOfl by a Taylor expansion of the coefhcients. 



III. A GENERALIZATION TO 1+1 
DIMENSIONS 

In this Section we generalize the OUP and construct 
in 1 + 1 dimensions a spatiotemporal random field <f{x, t) 
with reasonable properties. For fixed x the process 
should be the common temporal OUP described above 
and for fixed t a spatial OUP. It is reasonable to require 
translational invariance, analogous to temporal station- 
arity, of all averages and an exponential decay of the 
spatiotemporal autocorrelation 



Kf{x -x',t^ t') = {if{x, t)if{x',t')) 



^-\x-x'\/\-\t-t'\/7 



4Ar 



(13) 



where a = (TsCt, (Js a-nd A characterizing for fixed time 
the spatial process. 

We propose two independent schemes leading to the 
same result. First we employ a spatially discrete scheme 
to construct mo re g eneral spatiotemporal correlated pro- 
cesses given in '29\. Alternatively, we consider a linear 
stochastic PDE different from Eq. (^) driven by additive 
Gaussian spatiotemporal white noise and show that it's 
stationary solutions are Gaussian distributed and have 
the desired properties. The analytic results are corrobo- 
rated by numerical data. 



A. Recursive Generation 

We consider the field if{x,t) on equidistant lattice sites 
i, i = 0, . . . , N , adopting the notation ipi{t) = ip{Ax ■ 
i,t). In the first step of construction we generate spatially 
independent OUPs i]i{t) using a standard algorithm, e.g. 
|23.l27l| . with autocorrelation 



KUt-t') = {rn{t)vAt')) 



i_e"l«-«'l/Tx. . (lA) 



2 



We then, as proposed in [2g to construct a more gen- 
eral spatiotemporal correlated noise, superpose these pro- 
cesses 



(15) 



fe=0 



Since this expression is linear in rj^ also the (pi are Gaus- 
sian distributed with zero mean. Requiring that the spa- 
tial autocorrelation is the discrete version of Eq. for 
equal times, 



/c=0 



determines the coefficients ajk- It is easy to check that 



ajk 



, ^ ^-(j-k)Ax/X f _ ^-2Ax/\ 



(17) 



where j — Q, . . . ,N and k — 0, . . . , j. Using these coeffi- 
cients we write 



-jAx/\ 



Vo(t) 



+ ... + ^/T~P^^r^j{t). (18) 
With the corresponding formula for Lpj^i{t) we find 



where (i) are the spatially independent random num- 
bers specified above, each of which being a temporal 
OUP. Obviously, this is the spatial analogue of the dis- 
crete integral algorithm (|10|l for the temporal process. 

In discrete notation both for space and time we in- 
sert ryj+i_(+i from Eq. ifTIUl . after replacing fTt/V2r 
ct/V 4rA, into Eq. (|19|l (written for ipj+i^t+i) and obtain 
finally 



_g-At/r-Ax/A^^.^ 



-2A 



^AVl 



-2At/T. 



=5 



(20) 



where {£,i^tS,j,t') — ^i.jSt,t'- For simulations this discrete 
integral algorithm is preferrable since by construction it 
is correct for any choice of Aa; and At. 

Expanding the coefficients in Eq. 1)20(1 for small Ax 
and At a first order discrete differential algorithm, the 
generalization of Eq. p2|l . is obtained. 



^j+i,t+i 



1 - 



Ax 



At 
1 

T 



At\ 



AxAt ~ 
— ^^crCi+i,t+i, 



(21) 



where a = a / \fAxAi. Writing Eq. (|21() in continuous 
notation. 



Ax 

(p{x + Ax,t + At) ^ - — ] fix, t + At) 



1 - 



At 



ip{x + Ax, t) 



, Ax\ / AA , , 

^xAt^, ^ , , 

+^^£,{x + Ax,t + At), 



(22) 



we have replaced (r^j.t — > £,ix,t) in complete analogy 
with the rescaling for the temporal OUP. The autocor- 
relation of the spatiotemporal Gaussian white noise is 
(C(x, t)^{x',t')) = a^S{x - x')S{t - t'). 



B. Continuous Approach 

1. Generalized Langevin equation 

An alternative approach starts from a generalized 
Langevin equation in 1+1 dimensions, the stochastic 
PDE 



l + r| + AA + Ar^)^(x,t).C(.,.). (23) 



which for A = reproduces the Langevin equation (0) 
generating the temporal OUP and for r = that of the 
spatial OUP. The spatiotemporal Gaussian white noise 
can be conceived as the product ^{x, t) = ^{x)(,{t), where 
f (x) and ^(t) denote independent spatial and temporal 
Gaussian white noise, respectively. 

It is interesting to note that Eq. (|23|l is a hyperbolic 
PDE whereas Eq. J^l is a parabolic one. Eq. H23I) has the 
two families of characteristics x — const and t ~ const, 
the latter one is the only family of characteristics of Eq. 
(P) . Correspondingly, the solution of Eq. iQJ reproduces 
in the limit A ^ the temporal OUP multiplied by d{x — 
x') but r — > results not in the spatial OUP, see below. 

The stochastic PDE (I23|l can be obtained from the 
continuous differential algorithm (|22|l by Taylor ex- 
pansion of if for small Ax and At and performing 
the limit Ax, At 0. Alternatively, it can be 

conceived as the product of the two Langevin equa- 
tions for a temporal OUP, Eq. and its spatial 
analogue. For this we denote the product of the 
temporal and the spatial OUPs by ip{x, t) and ob- 
serve that the differential operator on the l.h.s. of Eq. 

factorizes as {l + rd / dt + \d / dx + Xrd"^ / dxdt) = 
{l + Td/dt){l + \d/dx). 

Using a separation ansatz, a solution of Eq. (|23|l can 



3 



be written as 



f{x,t) = f{x)g{t), where 
A 



(24) 



fix) = /oe-(---")/^ + -e--/^ / dW{y)ey^\ (25) 

g(i) = goe-(*-*°^/" + -pe-*/^ / dVK(s)e^/", (26) 

and /(xo) = /o and 5(^0) = 9o denote boundary and ini- 
tial values. The initial and boundary processes ip{x,t()) 
and ip{xQ, t) are OUPs with correlation length A and cor- 
relation time r, respectively. Note the appearance of 
the extra factors A and 1/A in Eqs. H25I) and (|26l) re- 
spectively, compared with the process given by Eq. Q). 
< 1^1 < c» is an arbitrary constant which corresponds 
to the separation constant for a deterministic PDE. In 
the nonstationary case it weights the relative influence 
of the initial and boundary realizations. In the term of 
^(x, t) which survives in the stationary case, A cancels 
and naturally its value plays no role, see below. 

Exploiting that the spatial and the temporal Wiener 
processes W{y) and W{s) are independent and have zero 
mean we obtain the autocorrelation function 

X'^(x,t; x',t') 



-\x-x'\/X-\t-t'\/T 

4tA 



0; 



1 ai 



-(x+x'-2xa)/X~\t~t'\/T 



+ ( {9l)A' 



2\ 42 ^ ] p-\x-x'\/\-{t+t'~2to)/r 

2t^ 

2 ^ 



xe 



12 2, 

(x+x' -2xo)/\-{t+t' -2tQ)/T 



^ A2 2A 



I1. 
2t 



(27) 



The first term on the r.h.s. is just the desired sta- 
tionary and homogeneous autocorrelation, independent 
on the boundary and initial conditions, cf. Eq. The 
remaining terms disappear for a;o — > — 00 and to —00, 
respectively. A second possibility to make the nonsta- 
tionary and nonhomogeneous terms vanish is to chose /o 
and go as zero mean Gaussian distributed with variance 
such that 



</»> - fx-"'- 



and 



W - 2r^2- 



(28) 



In this case the process will be homogeneous and station- 
ary from the beginning. The variance of the process f is 
independent of {x,t), {ip'^{x,t)) = (/o)(go) = o'^/4Ar. 



2. Generalized Fokker-Planck Equation 



spatially discretized system the FPE is a matrix equa- 
tion. We will derive for this case the generalized FPE 
and its stationary solution. Discretizing Eq. H23|l using 
a first order Euler-scheme gives the system of ordinary 
differential equations. 



rCdipit) = -Cip{t)dt ■ 



/Ax 



dW{t), (29) 



where <p — ((pi, . . . , i/jat)'^ and W = [Wi, . . . .Wn)^ ■ 
The matrix C has the non-vanishing elements 

Ci,i = Co = 1 + c-i+i.i = ci = (30) 
Ax Ax 

Since det C 7^ we can multiply Eq. H29I) with C^^ and 
obtain 



Td(p{t) — —ip{t)dt 



/Ax 



C-idW(t). (31) 



Now we can treat the system as a multivariate OUP. It 
can be shown that the corresponding FPE is 



(32) 



where p = p{ip, t\ipQ, to) is the transition probability den- 
sity and J is the hereby defined probability current den- 
sity. We note that in our case C^i {C-^)^ = {C'^C)^\ 
A stationary solution of Eq. (|32|l . means J = const. 
For natural boundaries where the probability current 
vanishes we have 



^ijfjPs 



2t2Ax ^ dip/' 



From Eq. 1)33(1 we get 

1- Inps = V 



2rAx 



2 



= 0. 

(33) 

(34) 



Now the non- vanishing elements of C^C can be com- 
puted from Eq. (|30|) as 



(C^C)^^^ ^cl + d (C^C) 



CqCi. 



(35) 



As the right hand side of Eq. H34|) is a gradient (C"'"C is 
symmetric), the potential conditions are fulfilled and a 
simple integration gives 



The Fokker-Planck equation (FPE) corresponding to a 
stochastic PDE should be a functional equation. For the 



Ps{ip) = TVexp 



(36) 



4 



where M is the normahzation factor. C'^^C is an oscilla- 
tion matrix with the positive eigenvalues 



A, 



A 



A2 



7V+ 1 



(37) 



Thus the stationary solution can be normalized, M 



Jl^j^ cr (tttAxAj) and the stationary probability 

density is indeed the zero mean Gaussian distribution 



C. Comparison with Numerics 

We compare the analytically given autocorrelation 
with numerically generated data obtained with the dis- 
crete integral algorithm provided by Eq. H20|l . Fig. 
shows a good agreement for fixed temporal and fixed spa- 
tial argument, respectively, imposing initial and bound- 
ary conditions which ensure stationarity and homogene- 
ity as described above. 

We also determined the mean square deviation of the 
variance of averages over 10^ independent realizations 
which is governed by the x^-distribution. The variance 
was always found within a 80 % confidence interval. 



IV. THE APPROACH OF GSR IN 1+1 
DIMENSIONS 

The above proposed generalization of the OUP has 
by construction autocorrelations decaying exponentially 
both in space and time. This is in contrast to the spa- 
tiotemporal correlated noise proposed by GSR IT3| . Since 
the autocorrelation for l-fl dimensions in real space was 
not explicitly calculated in the previous literature we be- 
low consider this case. Again, we derive the result in a 
continuous approach and in a spatially discretized scheme 
and compare the analytical results with numerical data. 
The autocorrelation in real space for spatial dimension 
d > 2 is evaluated in a different context in . In recip- 
rocal space, the result is given for general d in H^. 
cf. also lia. 



A. Continuous Approach 

We start with the Fourier transform of Eq. Q in d = 1 
which reads 

r^(p(/c, t) - -c(fc)^(fc, t) + t), (38) 

where c{k) = 1-1- A^fc^ and £,{k,t) is the Fourier trans- 
formed white noise with autocorrelation 

K^{k, t, k', t') = 2Tra^S{k + k')5{t - t'). (39) 



1.0 : 



0.1 



(a) 



0.001 
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FIG. 1: Autocorrelation of the generalized OUP in l-fl di- 
mensions normalized by the variance cr^/4Ar. Comparison 
of analytical and numerical results for (a) fixed temporal 
(t = t' = 100) and (b) fixed spatial {x = x' ^ 100) argu- 
ment. The lines show the analytic results from Eq. 1271 . the 
symbols are the results of simulations (squares: X = t — 100, 
triangles: A — 50, t = 80). Stationarity was ensured impos- 
ing the corresponding initial and boundary processes, see text. 
Averages over 10^ realizations (N=1000, a = 1, Ax = A/100, 
At = r/100). 



Equation (|38|l defines an OUP for each k. It has the 
general solution 



^{k,t) 



- p-c(fe)t/ 



Xfc,0) 



1 /■* 

_g-c(fc)t/r / ds^(fc,s)e=W''/^ 



(40) 



Stationarity and homogeneity is ensured if the initial val- 
ues have the autocorrelation 

K'^(k,0;k',0) = ^^Sik + k'), (41) 
2t c[k) 

as for the Fourier transform of a spatial OUP with vari- 
ance (t^/4tA. The autocorrelation function in the sta- 
tionary and homogeneous case is 



K'^{k,t;k',t') = 



(5(fc + fc')e-"('=)l*-*'l/^(42) 



which is up to constant factors in accordance with [13 . 
[Tgj i . Inverse Fourier transform gives 



K'^{x- x',t-t') 



2r 27r 



^c(fc) 



^-c{k)\t-t'\/T-ikix-x') 



, (43) 



5 



where we introduced Kf{x - x',t- t') = Kf{x, x'; t, t'). 
To calculate the integral on the right hand side of Eq. 
(I43|l we introduce k — Xk and c{k) = 1 + fc^ = c{k) 
and note that K"^ depends only on p = {x — x')/X and 
s = \t — t'l/r. The derivative of Eq. H43I) with respect to 
s reduces to the Fourier transform of a Gaussian 



dK'^ip, s) 



1 



a' 1 
2tA 27r 

-s— p^/4s 



dke 



4tA ^/^Fs 

Integration with respect to s gives 



(44) 



4rA 2 



e''erf(,+ ^ 



= ''erf w 



(45) 



where erf (a;) = die * is the error function. In the 

limit s — > oo the autocorrelation should vanish, hence 



lim - — - 

o^oo 8tA 




(46) 



erf 



P 



erf ( \/s — 



-erf 



'So 



The limit sg ^ oo should be carefully taken. If we are 
interested in the limit A — > or in the asymptotics for 
large p the corresponding operation has to be done before 
So ^ oo. The hmit A of Eq. (gS)) leads to 



K'^{s, X — x') = -^^{x ~ 3;')e 



(47) 



as to be expected. Evaluating first the limit r ^ of Eq. 
(|^ results in 



K^it-t',p)^-Sit~t')il + \p\)e-\''\. 



(48) 



Independent on the order of the limits we obtain for both 
A and r — > the result for spatiotemporal Gaussian white 
noise K'^ = cr'^5{x~ x')5{t — t') which can be also directly 
infered from Eq. 

The asymptotics for p 3> 1 and s = const ^ and, 
alternatively, for s ':$> 1 and p = const, is obtained em- 



ploying erf(z) 
[8fij |. as 



±1 



for 



±oo, cf. e.g. 



(49) 



For s <C 1 and p = const 7^ one obtains from Eq. 
(I46|l after first doing sq 00 and employing again the 
asymptotics of erf(2:) 



TT S — p2/4s 



(50) 



where the second term on the right hand side vanishes 
for s — > 0. 

For p <C 1 and s — const 7^ expanding 
erf (yi ± p/2y^) and e^P one obtains from inde- 
pendent on the order of the limits, up to second order in 
P 



(p,s)-^{(l-rfVi)(2V)-;^e-^}. (51) 



The limit r leads to K'P^p.s) = crV(4A)(5(t - 
t') (1 — p^/2) in accordance with the expansion of Eq. 
(I48|l for small p. Using the asymptotics 1 — erf-^s ^ 
e^^(\ — l/2s)/ y/¥s for large s one obtains K'f ^ 
CT^/(4TA-\/7ri)e~*'(l — p^ j^s) which agrees with the ex- 
pansion of Eq. (|49(l for small p. 

The autocorrelation function should solve the equation 
obtained by multiplying Eq. with ip[x' ,t') and aver- 
aging, 

T^K^ix -x',t^ t') = - (1 - A^A) K^{x ~x',t~ t') 



dt 

+e(i' - t) 



l_ 

2tA -v/Trs 



(52) 



which is fulfilled by 14t)|) . In the limit A ^ the inho- 
mogenity reduces to that of Q multiplied by 5{x — x') 
as it should be. In the limit t — > the inhomogeneity of 
Eq. ^ becomes a"^ / {2\)6{t - t')e-\p\ ^ which can be 
also directly derived. 

Sancho et al. 22] claimed that the decay of correlations 
is exponentially dominated in both space and time. The 
above results show that this is generally not the case for 
d—\, see 331. 



B. Spatially Discretized Scheme 

Garcfa-Ojalvo et al. 0| calculated the autocorrelation 
of the GSR process for d = 2 in discrete space. Here we 
repeat the procedure in d = 1 to compare it with the 
continuous case. The spatially discretized version of Eq. 
(P) reads 

T^^j(t) = -vAt) + A2A(p,(0 + i,{t), (53) 
where the Euler discretization of the Laplacian is 



A(p,(i) 



1 



Ax' 



(^j+i(t)-2^,(t) + ^,_i(t)). (54) 
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In discrete space we have to rescale the white noise ac- 
cording to 



(55) 



Again, as in the continuous case, we Fourier transform, 
solve the decoupled equations and calculate the autocor- 
relation function. We define the discrete Fourier trans- 
form on the spatial lattice as 



N-l 



/N)m 



Hence the inverse Fourier transform is given by 
1 



N-l 



-i(27r/Af)pj 



NAx ^ 



(56) 



(57) 



Greek indices are used in Fourier space and latin indices 
in real space. The indices run from to A^— 1 in both real 
and Fourier space; due to periodic boundaries — /i has to 
be interpreted as N — fj,. Now we can Fourier transform 
Eq. Q 



d 



where 



(58) 



(59) 



The autocorrelation function of the Fourier transformed 
white noise is 

(e^WC^'l^')) = 'j''NAx5^,-^,5{t - t'). (60) 



As in continuous space, Eq. H58() defines an Ornstein- 
Uhlenbeck process with autocorrelation time t/c^ for 
each /i. The stationary autocorrelation can be computed 
in complete analogy to continous space as the inverse 
Fourier transform of 



NAx 



2t 



Jti,-fj.' 



-C^\t-t'\/T 



(61) 



Hence the stationary autocorrelation in discrete space is 



Af-l 



1 V- 1 -c 
At ^ — ' r 



-c,,\t-t'\/r-i(2TT/N)tJ^(]-]') 



2t NAx ^ c„ 

M=0 



(62) 

Since = cm-^i, the imaginary part of the sum vanishes. 
For \t — t'\/T 3> 1 the autocorrelation is dominated by 
the first term e~l*~* '/'^ in the sum (/i = 0). However, 
this is not so for |i — t'| ~ r, cf. Fig.|2l3. 

Observing A^Aa; = L being the system size, and 
identifying k — 27r /i/L we have in the limit Ax — > the 
correspondence c{k) — \ + \^k^. Hence Eq. H61|) 

corresponds to H42|) and Eq. H62|) to H43|) after doing the 
limit i — > oo in an appropriate way. 




0.5 



0.0 




FIG. 2; Autocorrelation (normalized by the variance) of the 
GSR process in 1-1-1 dimensions. Comparison of simulations 
in Fourier space with analytical results, (a) shows the spatial 
dependence after a transient period t = t' — 1000 (note the 
symmetry due to periodic boundary conditions), and (b) the 
temporal dependence for x = x' = 32. Analytic results from 
the continuous approach (solid line), Eq. 14611 . and from the 
discrete approach (triangles), Eq. 16211 . practically coincide 
with numerical data (squares) obtained with Eq. Q. The 
dashed line in (b) shows for comparison the autocorrelation of 
a temporal OUP with r = 20. Averages over 10^ realizations 
(iV = 64, J = 1, Ax = = 1, A = 10, r = 20). 



C. Comparison with Numerics 

The initial conditions for a stationary field in Fourier 
space were chosen as independent Gaussian random num- 
bers with variance a'^{N Ax)^ /tc^ for each The spa- 
tial autocorrelation was computed using the correlation 
theorem (cf., e.g., [s^]) valid for weak stationary ergodic 
processes 



^ [{9{xo)9{xo + x))] = T[g{x)]T[~g{x)], 



(63) 



where J^[g(x)] denotes the Fourier transform of g{x). The 
procedure is faster in numerical simulations and gives the 
same results as the real-space approach, moreover the 
inverse Fourier transform can be avoided if one is only 
interested in spatial correlations. 

Figure 121 compares numerical and analytical results for 
the GSR process in l-fl dimensions. 

Simulations in real space give results which coincide 
with those in Fourier space; we refrain from demonstrat- 
ing this here. A simulation in real space has the disadvan- 
tage that the maximal possible time step At is restricted 
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5 10 15 20 25 

t - t' 



FIG. 3: Autocorrelation (normalized by the variance) of the 
GSR process in 1+1 dimensions. Comparison of simulations 
in real space for x = x' = 128 with periodic boundary con- 
ditions (triangles) and without periodic boundary conditions 
(squares). Both practically coincide with the analytical result, 
Eq. 1621 . Averages over IC realizations after a transient pe- 
riod t' = 500 (iV = 256, 5 = 1, Ax = 1, At = 0.1, A = 2, 
T = 20). 



by At < rArrfMA^, otherwise the discrete Eq. H53(l looses 
stabiUty, of. |13l |. 

Another issue is the dependence on boundary condi- 
tions which is shown in Fig. |31 Since in Fourier space 
we always have periodic boundary conditions we work in 
real space, where the above stability condition enforces to 
use a smaller At. In real space simulations without peri- 
odic boundary conditions one has to impose a stochastic 
boundary process. We used a temporal OUP with r = 20 
which has a different autocorrelation than the temporal 



autocorrelation of the GSR process. Therefore we show 
data after a transient period apart from the boundaries 
of the system. Clearly, all procedures give numerical data 
matching very well with the analytical result. 

V. CONCLUSIONS 

We introduced in 1+1 dimensions a spatiotemporal 
stochastic process with an autocorrelation exponentially 
decaying both in space and time, thus being a generaliza- 
tion of the OUP. An analogous generalization to higher 
spatial dimensions, although formally possible, seems 
physically not meaningful: The autocorrelation function 
should not factorize in the spatial variables. 

The situation resembles to that of the checkerboard 
process in 1+1 dimensions js^ [s^ Is^, 113 driven by a 
velocity changing randomly the sign which is modeled by 
the simplest discrete process with exponentially decaying 
autocorrelation, the dichotomous Markovian process. 
The checkerboard process is intimately connected with 
the Dirac equation or the Klein-Gordon equation in 
d — 1. Also there, the generalization to higher spatial 
dimensions meets nontrivial difficulties |36l |37| . 
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